Analytical results for coupled map lattices with long-range interactions 



O 

o 

(N 
G 

(N 



Q 

U 
d 



(N 
> 

o 

oo 
O 
m 
O 



> 
X 



Celia Anteneodo 1 , Sandro E. de S. Pinto 2 , Antonio M. Batista 3 and Ricardo L. Viana 2 
1, Centra Brasileiro de Pesquisas Fisicas, Rua Dr. Xavier Sigaud 150, 22290-180, Rio de Janeiro, RJ, 
2. Departamento de Fisica, Universidade Federal do Parana, 81531-990, Curitiba, PR, 
3. Departamento de Matemdtica e Estatistica, Universidade Estadual de Ponta Grossa, Ponta Grossa, PR, Brazil. 

(February 8, 2008) 

We obtain exact analytical results for lattices of maps with couplings that decay with distance 
as r~ a . We analyze the effect of the coupling range on the system dynamics through the Lyapunov 
spectrum. For lattices whose elements are piecewise linear maps, we get an algebraic expression 
for the Lyapunov spectrum. When the local dynamics is given by a nonlinear map, the Lyapunov 
spectrum for a completely synchronized state is analytically obtained. The critical lines character- 
izing the synchronization transition are determined from the expression for the largest transversal 
Lyapunov exponent. In particular, it is shown that in the thermodynamical limit, such transition is 
only possible for sufficiently long-range interactions, namely, for a < a c < d, where d is the lattice 
dimension. 

PACS numbers: 05.45.Ra,05.45.-a,05.45.Xt 



Synchronization between coupled chaotic systems is 
one of the most intriguing nonlinear phenomena [1] . It 
has attracted much interest since two decades ago [2] as 
it appears in a wide range of real systems such as in ar- 
rays of Josephson junctions [3], oscillating chemical reac- 
tions [4] , physiological processes [5] , and has applications 
as in communications [6] and control theory [7]. There 
are many types of synchronized behavior [8], but we are 
particularly interested in the completely synchronized 
states (CSSs) of coupled map lattices (CMLs), where all 
maps present the same amplitude at all times. Com- 
plete synchronization is an example of non-equilibrium 
phase transition [9] , which may be related to actual crit- 
ical phenomena like the superconducting-normal transi- 
tion in Josephson junctions [10]. 

CMLs, which are dynamical systems with discrete 
space and time, and a continuous state variable, have 
been investigated as theoretical models of spatiotemporal 
phenomena in a variety of problems in condensed mat- 
ter physics, neuroscience and chemical physics [11]. The 
spatiotemporal behavior is governed by two simultane- 
ous mechanisms: the intrinsic nonlinear dynamics of each 
map, and diffusion due to the spatial coupling between 
maps; the dynamical pattern being the outcome of the 
competition between them. This applies, in particular, 
to the problem of synchronization of chaotic maps [12]. 
The effective coupling range is a crucial factor to de- 
termine whether or not chaotic maps mutually synchro- 
nize. Nearest-neighbor couplings (short range) do not fa- 
vor synchronization, since the coupling effect is typically 
too weak to overcome the intrinsic randomness of map 
dynamics [13]. On the other hand, long-range couplings 
tend to facilitate synchronization, as exemplified by the 
limiting case of global (mean-field) coupling [14]. Lat- 
tices of non-locally coupled maps appear in neural net- 
works with local production of information [15], models of 
physico-chemical reactions [16], assemblies of biological 
cells with oscillatory activity [17], and diffusion coupling 
in nucleation kinetics [18]. Beyond CMLs, systems with 
many degrees of freedom with long-range couplings are 
an interesting object of study because of their anomalies 



(appearing at the level of the macroscopic thermodynam- 
ical description as well as in the underlying microscopic 
dynamics), which still require deeper understanding [19]. 
Simple dynamical models, such as CMLs, may add new 
knowledge on non-equilibrium long range systems. How- 
ever, there is a lack of analytical results for CMLs with 
arbitrary range couplings. Exact analytical results are 
particularly crucial because the occurrence of phenom- 
ena such as shadowing breakdown [20] or spurious syn- 
chronization [21] set difficulties in numerical approaches 
due to the unavoidable finite precision of numerical sim- 
ulations. 

Here we examine a form of coupling whose intensity de- 
cays with the distance r between sites as l/r a , with a > 
[22]. It has also been considered in biological networks 
[23], in ferromagnetic spin models [24], many-particle 
conservative (Hamiltonian time evolution) classical sys- 
tems [25,26], large populations of limit cycle oscillators 
[27] and a generalization of the Kuramoto model [28], 
among other examples. Explicitly, we consider a chain of 
N coupled one-dimensional chaotic maps x i— > f(x) such 
that the coupling prescription is 
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where x n represents the state variable for the site i 
(i = 1,2, ...,7V) at time n, e > and a > are the 
coupling strength and effective range, respectively, and 
rj(a) = 2^2^ = i r ~ a j is the normalization factor, with 
N' = (N — l)/2 for odd N. In conservative systems 
[25,26], scaling by 77 plays an important role in making 
the systems pseudo-extensive. Here periodic boundary 
conditions x^ — Xn and random initial conditions 
are assumed. The coupling term is a weighted average 
of discretized spatial second derivatives, the normaliza- 
tion factors being the sum of the corresponding statistical 
weights. It is straightforward to prove that in the limits 
a = and a — > 00 Eq. (1) reduces to the global mean- 
field and the local Laplacian-type couplings, respectively. 
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We characterize the spatio-temporal synchronization 
dynamics by means of the Lyapunov spectrum (LS) of 
the lattice, that enables one to estimate, for instance, 
the Kolmogorov-Sinai entropy through the Pcsin formula 
[29] and the Lyapunov dimension, which gives an up- 
per bound on the effective number of degrees of freedom 
needed to characterize the system dynamics [30] . Besides 
characterizing a CSS, when it exists at all, we must in- 
vestigate its stability with respect to small perturbations. 
If the CSS turns out to be dynamically unstable, we are 
faced with two possibilities: either the CSS presents the 
so called bubbling attractor, and in this case the CSS 
only lasts for a finite time, or the CSS loses transversal 
stability through a blowout bifurcation [8] . 

In this work we will present exact analytical results 
for the CML (1). We will show that for a ID lattice 
of N coupled piecewise linear maps it is possible to ob- 
tain an exact analytical expression for the LS, the results 
shown in [14,29] being recovered in the limits a — and 
a — ► oo. When the maps x i— » f(x) arc nonlinear, we will 
show that analytical results are still possible for CSSs. 
By means of the algebraic formulas for the LS, one can 
find the synchronization regions in the e x a space, since 
the second largest Lyapunov exponent (belonging to the 
direction transversal to the SM) equal to zero indicates a 
transition to the synchronized state [16,31]. Finally, the 
results obtained for a chain of maps will be extended to 
<i-dimcnsional hypercubic lattices. 

In order to calculate the LS one has to consider the 
tangent dynamics. By differentiating the equations of 
the original maps (1), one obtains the evolution equa- 
tions for tangent vectors £ = (Sx^, 5x^ 2 \ . . . , 5x^) T , 
that in matrix form read £ n+ i = T„£„, with the Jaco- 
bian matrix T„ given by 



T„ = 



(1-0 + 



T)(C 



-B 



(2) 



where the matrices D„ and B arc defined, respectively, 
by Di k = f'{x^)8 Jk and B jk = l/r%{\ - S jk ) , being 
tjh — min/g^lj — k + IN\. Notice that the particular 
choice of the interaction law is embodied in the matrix 
B which is time independent. 

Once specified the initial conditions, the LS is ex- 
tracted from the evolution of the initial tangent vector 
£,o- in = %i£,o, where T n = T„_i . . . TiT is product 
of n Jacobian matrices calculated at successive points of 
the discrete trajectory. If Ai, . . . , Ajv are the eigenvalues 
of A = lim (T^T n )2^ (that are real and positive), the 

n— *oo 

Lyapunov exponents are obtained as [32] 



A fc =lnA fe , k = l,...,N. 



(3) 



We start by applying the expressions above to the 
piecewise linear maps x t— > f(x) = (3x (mod 1), with 
(3 > 1. In this case we have f'(x) = (3 = constant, there- 
fore D„ = /31jv, and T„ becomes 
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the rightmost identity defining the matrix B. Since the 
symmetric tangent map does not depend on time, it re- 



sults A = /3B. So, in order to obtain the LS, it is enough 
to diagonalize B. Because of its periodicity, B can be 
diagonalized in Fourier space [33], the eigenvalues being 



N' 



cos(2nkm/N) 



, N, 



m—l 



(5) 



where we considered odd N. Finally, from Eq. (3), tak- 
ing into account the special form of A, the LS is given 

by 



A fe = ln/3 + ln 



l-e + 
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This expression is consistent with previous numerical re- 
sults [34]. In the extreme cases a — > 00 and a = the 
known expressions [14,29] are recovered. 

Now we will consider lattices of nonlinear maps. An 
important case that can be tackled easily is the one where 
the maps are in the CSS. As it will become clear soon, this 
instance provides relevant information on the synchro- 
nization transition. In the CSS, the dynamical variables 
of all maps coincide, i.e, Xn^ — x^ = ... ,x^ = Xn\ 
at each time step n. The LS for the CCS when a = has 
already been found by Kaneko [14]. Now, for arbitrary 

a, we have D„ = f'{xn^)^-N, thus, T„ = f'{xn^)~B and 
TjTn = {I\%l[f{x { * ] )] 2 )B 2n . Therefore, for the CSS, 
following Eq. (3), one arrives at 

\l = lim J2ln\f'(x\* ) )\+ln\l-e+-^b k \, (7) 

n^oo n £ -~' 77 a 



i=l 



where b k are the eigenvalues of B defined in (5). Assum- 
ing ergodicity, the time-average in (7) can be substituted 
by an average over the single-map attractor. In this way 
one gets 



\* k = \ v + In 



!-£ + 
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bk 



(8) 



where \jj — (In is the Lyapunov exponent of 

an uncoupled map. This expression is general: it ap- 
plies to any lattice of nonlinear ID maps coupled with 
the scheme here considered, the parameters that define 
the particular uncoupled map affecting only Xjj. For in- 
stance, for the logistic map x ^ f{x) = ax(l — x), with 
a = 4 and x G [0,1], \ v = (In |4(1 - 2x^)\) = In 2 [35] 
and the contribution of the power-law coupling is always 
ln|l — e + -^)b k \. Notice that the LS in CSSs has the 
same structure as the LS obtained for piecewise linear 
maps [Eqs. (5)-(6)]. 

The synchronization transition can be characterized 
by a complex order parameter [4] defined, for time n, 
as R n = \j; X^jLi e 2 ™"' |. A time-averaged amplitude 
R is computed over an interval large enough to warrant 
that the lattice has attained the asymptotic state. In the 

CSS, one has R = l. On the opposite case of completely 

(?') 

non-synchronized maps, the site state variables 
so uncorrelatcd that R ~ 0. 
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A diagnostic of synchronization can also be extracted 
from the LS. It can be easily verified that, for arbitrary 
a, the CSS lies along the direction of the eigenvector as- 
sociated to the largest exponent. This was previously 
observed by Kaneko for the particular case a = [14]. 
Therefore, the CSS will be transversally stable if the 
(N — 1) remaining exponents are non positive, that is, 
< (where the tilde stands for ordered exponents). 
The second largest exponent, corresponds in Eq. (8), 
to k = 1 (or k = N — 1, due to degeneracy) if the ar- 
gument of the modulus in Eq. (8) is positive and to 
k = (N — l)/2 (or k = (N + l)/2) otherwise. One ob- 
tains 



(1 -c- Xu )[ 1 



e' c = (l + e- x ") 1 



N' 

— Y 



cos{2irm/N)\ 1 



(9) 
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where e c (e^, ) are the coupling strengths below (above) 
which the SM ceases to be transversely stable, such that 
synchronized states are not typically observed. It is note- 
worthy that Eq. (9) is quite general for the coupling 
scheme here considered: the parameters that define the 
particular uncoupled nonlinear map (embedded in Xjj) 
just participate through the first factor. 




a 

FIG. 1. Synchronization in parameter plane e x a for a ID 
lattice of N coupled logistic maps x 1— > f(x) = 4x(l — x) and 
e < 1. The critical line e c (a) was determined analytically 
from Eq. (9) for different values of N (full lines) and numer- 
ically from the conditions E / 1 (open symbols) and \% = 
(full symbols) for N = 21 and 51. 

Fig. 1 presents the lower critical line in the parameter 
space e x a for N coupled logistic maps x 1— ► f(x) = 
4x(l — x). The critical line was obtained analytically 
from Eq. (9) and numerically by means of two proce- 
dures: either when R ^ 1 is numerically detected or by 
the condition of nullity of the second largest Lyapunov 
exponent in the CSS. The numerical results are in good 
agreement with the analytical prediction. 



As can be observed in Fig. 1, the critical frontier de- 
pends on the system size N . In the limit N — » 00 we 
obtain 



1 - e~ AL ' 
1 - C{a) ' 

where C(a) corresponds to 



(10) 



C(a) 



lim 



r N' 

E 



cos(2irm/N) 



N' 

E 



1 

m° 



(11) 



This limit is equal to unity for a > 1, so that Eq. (10) 
furnishes a divergent result. For a outside the domain of 
convergence of the series, one gets 



1-a r 
a ) = — i — / 

7T Jo 



cos(x) 



dx. 



(12) 



Moreover, one has e' c ^ = 1 + c~ Xu , if a < 1. Then, 
in the limit N — > 00, synchronization is only possible for 
sufficiently long-range interactions (see Fig. 2), namely, 
for a < a c < 1. The critical value here obtained for 
ID CMLs is different from the one reported for other 
ID systems with similar power-law interactions, such as 
ferromagnetic spin models [24] or many-particle classical 
Hamiltonian systems [25]; in such cases, the critical value 
for the existence of an order/disorder transition is a c = 2. 
In a generalization of the Kuramoto model, the value of 
a c is controversial: While a c = 2 was first reported [28], 
recent analytical considerations point to a c = 1 [36] . Al- 
though the generalized Kuramoto model is a continuous 
time dynamical system, our analytical result for the up- 
per bound of the critical value suggests that the latter 
value is the correct one. 

The expressions we have derived for the LS of chains 
of maps can be straightforwardly generalized for hyper- 
cubic lattices of arbitrary dimension d. In fact, for the 
general (i-dimensional case, it is straightforward to show 
that the eigenvalues of the matrix B become 



b k = 2 d ]T 



cos{2irrk-m/Nd) 



k=l,...,N, (13) 



where is the position vector of site k, m = 
(mi, . . .,m d ), with < m t < N', and N' = (N* - l)/2. 
The normalization factor reads 



r)(a,d) 



E- 



N a/d _ ^ 

1 - a/d ' 



(14) 



Then, following the same lines as before but now for 
the sake of generalizing Eq. (10), it is easy to see that 
e c ,oo will diverge if a/d > 1, which leads to a c < d. 

Additionally, our results for the LS could be extended 
to the more general class of coupling schemes where the 
dependence of the coupling strength on the inter-map dis- 
tance is not necessarily of the power-law type. In these 
cases, one should feed Eqs. (6) and (8) with the eigen- 
values of the appropriate matrix B, that contains the 
particular dependence of the interaction strength on dis- 
tance. 
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FIG. 2. Critical coupling strengths for synchronization vs 
range parameter, for different values of N indicated in the 
figure. The critical lines were analytically obtained from Eq. 
(9) for a ID lattice with power-law couplings and Xu = In 2. 
Synchronization occurs for e c < £ < e' c . 

In conclusion, we have presented analytical expressions 
for the LS of CMLs with an interaction which decays 
with the lattice distance as a power law, for two cases: 
(i) piecewise linear coupled maps; and (ii) the CSS of 
lattices of one-dimensional maps. Our results enable us 
to predict the critical values for synchronization in the 
coupling parameter plane (strength versus range), and 
also may be used to obtain related quantities of inter- 
est, like KS-entropies and Lyapunov dimensions. Such 
exact analytical results are crucial in order to avoid dif- 
ficulties present in numerical approaches, such as shad- 
owing breakdown, due to unavoidable finite precision of 
numerical simulations. In addition, we have shown that, 
in the thermodynamical limit, the critical range for syn- 
chronization is equal to the lattice dimension. Many of 
our results could be extended to lattices of continuous- 
time oscillators, and hence have an even wider range of 
applicability. 
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